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Evidence of strong electron-phonon interaction in superconducting MgB2 from 

electron tunneling. 

A. I. D'yachenko,[] V. Yu. Tarenkov,[] A. V. Abal'oshev, and S. J. Lewandowski 
Instytut Fizyki Polskiej Akademii Nauk, Al. Lotnikdw 32, 02-668 Warszawa, Poland 

We report tunneling measurements of the electron-phonon (e-ph) interaction in superconducting 
MgB2 using the MgB2-I-Nb junctions, where I stands for insulator. The phonon structure in tun- 
neling density of states in MgB2 clearly indicates strong e-ph coupling for the E2 9 in-plane boron 
phonons in a narrow range around 60 meV. The Eliashberg spectral function o? (w)F (u>) recon- 
structed from the tunneling data, exhibits significant additional contribution into e-ph interaction 
from other vibrations such as acoustic (~ 38 meV) and optical (~ 90 meV) bands. Our results are 
in reasonable agreement with neutron scattering experiments, and also to some data of Raman and 
infrared spectroscopy. 

PACS numbers: 74.25.Kc, 74.50.+r, 74.70.-b 



The recent discovery of superconductivity near 40K 
in MgB2, stimulated a renewed interest in the pair- 
ing mechanism of cuprate superconductors. It is gen- 
erally believed that understanding superconductivity in 
the simple binary MgB2 compound should be much eas- 
ier than in the high-T c superconductors. Theoretical ab 
initio calculations show that dominant contribution to 
Eliashberg function in MgB2 arises from the Raman- 
active E2 S phonons near the Brillouin zone center at 
~ 75 meV jj], ||, ||, 0, 0L This conclusion is supported 
by Raman studies @> which indicate a mode at 

~ 600 — 620 cm -1 (74 — 77 meV) with a very large width 
~ 200 cm -1 (25 meV), usually related to anharmonic ef- 
fects. However, there has been no consensus about the 
energy and role of the E2 3 mode, because this mode has 
not been identified unambiguously by inelastic neutron 
scattering or in Raman experiments ^ |], |l(| [H], |l2| ]. 

Tunneling spectroscopy can provide a direct mea- 
sure of the Eliashberg electron-phonon spectral function 
a 2 (oj)F(uj), which can reveal the importance of differ- 
ent phonon modes in the superconductivity of MgB2. 
Here F(u;) is phonon density of states and a(ui) is an 
effective electron-phonon coupling function for phonons 
of energy hw. Tunneling conductivity a = dl /dV gives 
direct information about the density of states N(lo) oc 
Re{w/[w 2 - A(w) 2 ] 1 / 2 } related to quasi-particle excita- 
tions in the superconductor, where A(u/) is a complex 
and energy dependent gap function Jl5| . a 2 (u})F(ui) can 
be obtained by the inversion of the Eliashberg equations 
for A(a->) 0. The method is based on the determina- 
tion of A(cj), therefore it non-sensitive to the presence of 
non-superconducting impurity phases in MgB2 . The lat- 
ter can give comparable features in Raman and infrared 
spectra at energy range typical for the superconducting 
MgB 2 @0. ' 

We report here high resolution superconducting tun- 
neling spectroscopy measurements of the phonon struc- 
ture of MgB2. We used tunnel point-contact junctions 
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FIG. 1: Experimental conductance of MgB2-I-Nb tunnel junc- 
tion in superconducting state (solid line) and the correspond- 
ing approximation of the normal state data (dashed line). 



Nb-I-MgB2, where I is a natural tunnel barrier. We found 
that the second derivative of tunnel current d 2 I/dV 2 ex- 
hibits features, which coincide with peaks in the phonon 
density of states (DOS) F(u) obtained from neutron scat- 
tering Raman g2j and infrared |0| [b§ spec- 
troscopy. Numerical treatment of the experimental data 
indicates that the observed phonon structures between 
55 — 65 meV and near 90 meV are very strongly coupled 
to the electrons. 

The investigated samples were prepared by compact- 
ing commercially available MgB 2 powder (Alfa Aesar, 
purity 98%) at high pressure (P > 5 GPa) into thin 
rectangular bars of about 1 x 0.6 x 0.08 mm 3 . Details 
of the sample preparation are described elsewhere 
The samples revealed a superconductive transition tem- 
perature (onset) at T c = 38 K and a high critical cur- 
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rent density J c > 10 5 A/cm 2 at T = 4.2 K, considerably 
larger than the current densities employed in the tunnel- 
ing measurements. Therefore, current depairing had no 
effect on the phonon-mediated region of tunneling data. 
The point-contact junctions were obtained by pressing 
Nb wire (0 ~ 50^*m) into the surface of the sample. The 
tunnel characteristics were measured using a standard 
four probe lock-in technique. 

In a number of Andreev S-N-S and tunnel S-I-S ex- 
periments at T = 4.2 K we observed two gaps A, the 
smaller gap value of about 3 meV, and a larger gap close 
to 7 meV, both results in agreement other measurements 
p2| . In tunnel S-I-S junctions we observed usually only 
the smaller gap with s-wave symmetry. Let us observe 
that the smaller ~ 3 meV and the larger ~ 7 meV gaps 
are associated with the 3d parts and quasi-2D sheets of 
the Fermi surface, respectively [Q. The experimental de- 
tection of the single small gap is possible in the case of 
tunneling perpendicular to the honeycomb born planes. 
Due to the strong directionality of the tunnelling exper- 
iments, the probability of normal injection is maximal 
and electronic properties are mainly probed along the 
direction perpendicular to the junction surface, i.e. in 
our case normal to the born planes. 

At high bias voltages (> 20 meV) the dI/dV(V) data 
reveal a structure (Fig.l), which in the case of conven- 
tional superconductors is characteristic of phonon effects 
p3[ . To investigate this structure in more detail, we used 
the second derivative of the tunneling current d 2 I/dV 2 , 
obtained by standard harmonic detection, as well as by 
straightforward numerical differentiation of the tunnel- 
ing conductivity dl/dV. Both methods give analogous 
results. In these spectra, we observed sometimes peaks 
at ~ 17 meV and ~ 28 meV, which are not intrinsic 
to clean MgB2. The structure near 28 meV most likely 
arises from a thin pure Mg layer in proximity with the 
superconducting MgB2 (analogous structures have been 
observed in Ag-I-Mg/Nb junctions f24[). For this rea- 
son, we excluded these singularities from further consid- 
erations, especially as they do not affect the subsequent 
treatment of the raw data. 

The plot d 2 I/dV 2 (V), shown in Fig. 2, at eV > 30 meV 
directly reflects the electron-phonon spectral function 
a 2 (w)F(w) of MgB 2 ; in particular the dips in d 2 I/dV 2 
correspond to peaks in a 2 (uj)F(uj). Peaks in F(lu) should 
be reproduced in a 2 (u;)F(w) at approximately the same 
energy, even if the coupling factor a 2 {u>) is strongly en- 
ergy dependent. The dl/dV and d 2 I/dV 2 data were used 
to obtain the normalized tunneling conductance a /aw, 
where un is a normal-state (background) conductance, 
and the tunneling density of states N(oj); the relevant 
numerical method is given in [23|j . 

The normalized density of states n(u) = 
A r (w)/7VBCs(w) — 1, proportional to the smoothed 
dl/dV curve, is shown in Fig. 3. Here -/Vbcs( w ) is 
the broadened BCS density of states (Dynes function). 
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FIG. 2: Second derivative d 2 I/dV 2 plots of a MgB 2 -I-Nb 
junction. The minima in the curve correspond to peaks in 
Eliashberg function a 2 (lo)F(uj). Voltage origin is at Am s b 2 + 
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FIG. 3: The tunneling density of states of MgE>2 n(uj), nor- 
malized to the BCS density of states. Experimental data are 
plotted as a solid line; dashed line shows n(u>) calculated from 
Eliashberg function a 2 (w)F(o;). 

The dl/dV curve displays a well-known behavior due 
to the phonon mode reflection in tunneling S-I-S con- 
ductance |l5|, [l(|, namely, sharp drops of dl/dV near 
the phonon emission thresholds, to which correspond 
dips in d 2 //dV^ 2 . We can rule out very weak inelastic 
process of phonon-assisted tunneling, because such 
process increases the tunneling conductance near the 
peaks in F(w) and gives corresponding peaks in d 2 I/dV 2 
p3[ . The common feature of all measured d 2 I/dV 2 
curves is the reproducibility of dip positions at ~ 40, 
55 -7- 65 and 90 -j- 92 meV, measured from the gap sum 
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Am 9 b 2 + ^Nb (Fig. 2). There is a close correlation in 
the peak positions in our d 2 I/dV 2 data and thouse in 
the generalized phonon density of states (GPDOS) as 
measured by inelastic neutron scattering data |ll| [l9| . 
The first tunneling peak near 40 meV arises from 
acoustic modes of MgBa. The well- resolved tunneling 
structures at ~ 60 meV and ~ 90 meV are associated 
with optical vibration of the born atoms. As it is well 
known, in the vicinity of peaks in a 2 (w)F(o;), the gap 
function A(u>) has a strong energy dependence and 
large Im{A(w)} [jl5|. At fku 3> |A| tunneling density 
N(u) ~ 1 + (l/2)Re{A 2 (w)/w 2 }, so the phonon struc- 
ture in N(oj) is weighted by the factor (SA/lu) 2 . In the 
strong coupling limit (the McMillan's electron-phonon 
coupling constant A of the order of 1) variation of A(uj), 
SA ~ A, and the deviation in N(u>) is of the order 
1 — 2% (Fig. 3). Even without detailed calculations one 
can see that since the amplitude of tunneling structure 
at 40 meV is roughly one half of that at 60 meV (c.f. 
Fig. 2) and the optical phonon energy is 1.5 times higher, 
the actual electronic coupling to the phonons at 60 meV 
may be as much as (1.5) 2 ~ 3 times greater than that 
for acoustic phonons at ~ 40 meV. On the other hand, 
it is evident, that the optical modes at ~ 90 meV give 
approximately the same contribution than that of modes 
at 60 meV. Such conclusion, obtained directly from the 
tunneling data, is independent of the anisotropy of the 
energy gap A(k) and agrees well with similar estimations 
based upon Raman and infrared spectra measurements 

An exact determination of a 2 (uj)F(u>) requires full in- 
version of the Eliashberg equations. We used normalized 
conductivity o~/<jn for self-consistent calculation of the 
spectral function a 2 (u>)F(uj) with a standard inversion 
program p3], which allows to take into account the al- 
ready mentioned possible presence of a thin Mg layer on 
the surface of MgB2, in which reduced superconductivity 
is induced by proximity effect. The gap edge region below 
26 meV was cut off to avoid a divergence of DOS. The 
obtained "experimental" function a 2 (oj)F(o;) is shown in 
Fig. 4 together with theoretical ab initio phonon density 
of states F t (uj) combined with neutron scattering jL9]. 
In isotropic limit for A = 0.9 and /i* = 0.1 Eliashberg 
function yields T c « 38 K, and |A| = 6.5 meV. Neutron 
GPDOS j£J[ H differs in its detailed shape from F t (w), 
but peaks in F t (w), GPDOS and our a 2 (w)F(w) all occur 
at the same energies. They all show pronounced minima 
at 40 -j- 50 meV and maxima at 37, 62 and 90 -j- 92 meV. 
The peak in a 2 (w)F(w), responsible for acoustic bands 
(~ 38 meV), is strongly suppressed, but does not disap- 
pear completely, as in the theoretical calculations |l|, Q. 
We do not see also the prevailing amplitude of one mode 
at 60 — 70 meV, as predicted for the theoretical function 
a 2 (cj)F t (cj) |l} |[ ||, although our experimental data indi- 
cate two bands of optical modes corresponding primarily 
to the born motions at ~ 60 and ~ 90 meV. For compar- 
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FIG. 4: Comparison of the phonon density of states in MgB2 
obtained by different methods, (a) The function a 2 (uj)F(lo) of 
MgB2 , determined from tunneling measurements in this work 
(open circles). The continues curve is the phonon density of 
states Ft(o») calculated from neutron data p9|| . The dotted 
curve shows Raman spectra of MgB2 at room temperature 
[[[ill , for the sample obtained under hig h pressure, (b) The 
calculated phonon dispersion curves 



ison, in Fig. 4(b) the calculated phonon dispersion curves 
are shown along the high-symmetry directions of the Bril- 
louin zone fll9"| . As seen, the Ei„ and A2 U infrared active 
modes give noticeable contribution to the a 2 (oj)F(uj) at 
40 and 50 meV, respectively. They have been also ob- 
served in infrared absorption experiments at 333 cm -1 
(41 meV) and 387 cm" 1 (48 meV)@. Optical studies 
in the c-axis oriented superconducting MgB2 films show 
two strong phonon peaks at 380 cm" 1 (48 meV) and 480 
cm" 1 (60 meV) jl3|. The peak in the "experimental" 
a 2 (u})F(uj) at 60 -j- 65 meV arises mainly from the E2 S 
phonon modes with wave vector along the T— A direction 
[Fig. 4(b)]. This phonon mode is Raman-active. In most 
investigations the Raman scattering in MgB2 revealed a 
broad peak near 620 cm -1 (78 meV) but, Meletov et al 
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|L2| have observed a peak near ~ 590 cm -1 (73 meV). 
They investigated ceramic MgB2 samples using a micro- [1 
Raman system, which rendered possible the identification 
of small crystalline grains of MgB 2 , whose Raman spec- 1- 
tra differ drastically from those of MgO or metallic Mg 
inclusions. After homogenisation of ceramic samples in- 
duced by high pressure (> 5 GPa) the E 2g phonon mode ; 
was observed at 500 — 550 cm -1 (62 — 68 meV) [the dot- 
ted curve at Fig. 4(a)], which is close to our a 2 (w)F(cj) [5 
result. For the E 2g mode, there are also large discrep- 
ancies between the results of various calculations. For " 
example, the harmonic energy w#(E 2g ) = 60 meV p9), 
but anharmonicity leads to an increase of wh(E2 S ) up 
to 75 meV flj ^9). Approximately the latter value (75 jy 
meV) was reported in the Raman study fij, where E 2g 
mo de had very large (~ 25 meV) width. On the other [8 



hanji, our experimental a (u)F(uj) coincides best with 



the harmonic value uoh (E 2g ) ~ 60 meV |19(] (Fig. 4) and 
is reasonably narrow. The extremely broad E 2g line at 
~ 620 cm -1 (78 meV) can include a sum of peaks from 
impurity phases, like Mg, MgO (phonon mode of MgO 
at ~ 80 meV has been observed in Ref. f25|), B 2 03, 
and second-order Raman signal from the acoustic ~ 300 



H9 

[10 
[11 



cm 1 (37 meV) phonon [ 10 . On the other hand, the 



infrared absorption spectrum of MgB 2 is characterized 
by a broad band centered at 485 cm -1 (60 meV) [ O , 
which can be associated with the peak in phonon density 
of states F(ui) at 60 meV, in agreement with our results 
(Fig.4). 

In summary, we have obtained high resolution tunnel- 
ing spectra of MgB 2 with T c = 38 K. The tunneling den- 
sity of states deviated from the standard Dynes function 
above 30 meV in a way, which is characteristic of phonon 
effects. The second derivative curve displays dips at en- 
ergies, which correspond to peaks in the phonon density 
of states measured by neutron scattering an also to some 
data of Raman and infrared spectroscopy. The spectral 
function a 2 (w)F(w) has been reconstructed from tunnel- 
ing data. The relatively big structure in tunneling spec- 
tra confirms strong coupling of electrons to optical E 2g 
phonon modes at ~ 60 meV. However, other vibrations 
such as acoustic (~ 38 meV) and optical (~ 90 meV) 
bands are also important and give noticeable peaks in 
the tunneling spectra. Therefore, the predominance of 
the E 2g phonon mode near the T point in the electron- 
phonon mechanism of superconductivity in MgB 2 may 
be questioned. 
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